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We have studied the phase diagram and entanglement of the one dimensional Ising model with 
Dzyaloshinskii-Moriya (DM) interaction. We have applied the quantum renormalization group 
(QRG) approach to get the stable fixed points, critical point and the scaling of coupling con- 
stants. This model has two phases, antiferromagnetic and saturated chiral ones. We have shown 
that the staggered magnetization is the order parameter of the system and DM interaction produces 
the chiral order in both phases. We have also implemented the exact diagonalization (Lanczos) 
method to calculate the static structure factors. The divergence of structure factor at the ordering 
momentum as the size of systems goes to infinity defines the critical point of the model. Moreover, 
we have analyzed the relevance of the entanglement in the model which allows us to shed insight 
on how the critical point is touched as the size of the system becomes large. Nonanalytic behavior 
of entanglement and finite size scaling have been analyzed which is tightly connected to the critical 
properties of the model. It is also suggested that a spin-fluid phase has a chiral order in terms of 
new spin operators which are defined by a nonlocal transformation. 

PACS numbers: 75.10.Pq, 73.43.Nq, 03.67.Mn, 64.60.ae 



I. INTRODUCTION 

At zero temperature, the properties of quantum many- 
body system is dictated by the structure of its ground 
state. The degree of complexity of this structure is dif- 
ferent for various systems. It ranges from exception- 
ally simple case (when a strong magnetic field aligns all 
the spins of a ferromagnet along the field direction, pro- 
ducing a product or unentangled state) to more intri- 
cate situation where entanglement pervades the ground 
state of system. Thus, entanglement appears naturally 
in low temperature quantum many body systems, and 
it is at the core of relevant quantum phenomena, such 
as superconductivity 1 , quantum Hall effect^,, and other 
quantum phase transitions 3 -. Quantum phase transitions 
have been one of the most interesting topics of strongly 
correlated systems during the last decade. It is basically 
a phase transition at zero temperature where the quan- 
tum fluctuations play the dominant rolei Suppression of 
the thermal fluctuations at zero temperature introduces 
the ground state as the representative of the system. The 
properties of the ground state may be changed drastically 
shown as a non-analytic behavior of a physical quantity 
by reaching the quantum critical point (QCP). This can 
be obtained by tunning a parameter in the Hamiltonian, 
such as the magnetic field or the amount of disorder. The 
ground state of a typical quantum many body system 
consists of a superposition of a huge number of product 
states. Understanding this structure is equivalent to es- 
tablishing how subsystems are interrelated, which in turn 
is what determines many of the relevant properties of the 
system. In this sense, the study of entanglement offers an 
attractive theoretical framework from which one may be 



able to go beyond customary approaches to the physics 
of quantum collective phenomena-. 

Recently some novel magnetic properties were dis- 
covered in a variety of quasi-one dimensional materials 
that are known to belong to the class of Dzyaloshinskii- 
Moriya (DM) (D .(Si x Sj)) magnet to explain helical 
magnetic structures. The relevance of antisymmetric su- 
perexchange interactions in spin Hamiltonians which de- 
scribe quantum antiferromagnetic (AF) systems was in- 
troduced phenomcnologically by Dzyaloshinski: 6 . Moriya 
showed later, that such interactions arises naturally in 
the perturbation theory due to the spin-orbit coupling 
in magnetic systems with low symmetry 7 -. Some AF 
systems are expected to be described by DM interac- 
tion, such as Cu(C 6 D 5 COO) 2 3D 2 0&2., Yb^As^^, 
BaCu 2 Si 2 7 1? , a - Fe 2 3 , LaMnO^- and K 2 V 3 Og£; 
exhibit unusual and interesting magnetic properties due 
to quantum fluctuations and/or in the presence of an 
applied magnetic fiel d 14 : 16 ' 17 . La 2 CuOi also belongs 
to the class of DM antiferromagnets, which is a parent 
compound of high-temperature superconductors^. This 
has stimulated extensive investigations of the properties 
which are created from DM interaction. This interaction 
is however, rather difficult to handle analytically, which 
makes the interpretation of experimental data to be hard. 
In addition, more knowledge in this respect expand our 
understanding of many interesting quantum phenomena 
of low-dimensional magnetic materials. 

Recent discovery of an unusual strong coupling be- 
tween the ferroelectric (FE) and magnetic order param- 
eters has also revived the interest in the magnetoelec- 
tric effecti 1 ^. Due to the possibility of easily controlling 
the electrical properties using magnetic field, search of 
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eigenstates of the blocks the full Hamiltonian is projected 
onto these eigenstates which gives the effective (renor- 
malized) Hamiltonian. 

The Hamiltonian of Ising model with DM interaction 
in the z direction on a periodic chain of N sites is 



FIG. 1: (color online)The decomposition of chain into three 
site blocks Hamiltonian (H B ) and inter- block Hamiltonian 
(H BB ). 



compounds, in which the magnetic order is incommen- 
surate with lattice period, is of particular interest for 
future applications32i21. Generally, certain types of mag- 
netic order can lower the symmetry of the system to 
that of the polar groups, which allow for ferroelectric- 
ity. According to the recent experimental results, he- 
lical magnetic structure are the most likely candidates 
to host ferroelectericit y 22 ' 23 i 24 . It has been shown that 
the DM interaction induces a FE lattice displacement 
and helps to stabilize helical magnetic structures at low 
temperature 2 ^. 

In the present paper, we have studied a one dimen- 
sional AF Ising model with DM interaction using the 
quantum renormalization group (QRG) and numerical 
exact diagonalization methods. In the next section the 
QRG approach will be explained and the renormaliza- 
tion of coupling constants are obtained. In section III, 
we will obtain the phase diagram, fixed points, critical 
points and the staggered magnetization as the order pa- 
rameter of the model. We will also introduce the chiral 
order as an ordering which is created by DM interaction. 
The exponent which shows the divergence of correlation 
function close to the critical point (u) , the dynamical ex- 
ponent (z) and the exponent which shows the vanishing 
of staggered magnetization near the critical point will 
be also calculated. We then present the numerical exact 
diagonalization results on finite sizes of N = 12, 16, 20 
and 24. In section IV we will calculate the renormal- 
ization of entanglemen t 26 ' 27 for this model and we will 
show that it has a scaling behavior near the QCP which 
is directly related to critical properties of the model. 



II. QUANTUM RENORMALIZATION GROUP 

The main idea of the RG method is the mode elim- 
ination or the thinning of the degrees of freedom fol- 
lowed by an iteration which reduces the number of vari- 
ables step by step until a more managable situation is 
reached. We have implemented the Kadanoff 's block ap- 
proach to do this purpose, because it is well suited to 
perform analytical calculations in the lattice models and 
they are conceptually easy to be extended to the higher 
dimension a 28 ' 29 ' 30 ' 31 . In the Kadanoff 's method, the lat- 
tice is divided into blocks in which the Hamiltonian is 
exactly diagonalized. By selecting a number of low-lying 
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The effective Hamiltonian (in the first order renormal- 
ization group prescription) is 
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We have considered a three site block procedure defined 
in Fig.©. The block Hamiltonian (H B = J2 h f) of tne 
three sites and its eigenstates and eigenvalues are given 
in Appendix A. The three site block Hamiltonian has 
four doubly degenerate eigenvalues (see appendix A). Pq 
is the projection operator of the ground state subspace 
which is defined by (Pq = | fr)(i/>o| + | -H-X^ol)) where 
\ipo) and \ijj' ) are the doubly degenerate ground states, 
| fr) and | -IJ.) are the renamed base kets in the effective 
Hilbert space. We have kept two states (\ipo) and \ip' )) 
for each block to define the effective (new) sites. Thus, 
the effective site can be considered as a spin 1/2. The 
effective Hamiltonian is not exactly similar to the initial 
one, i.e, the sign of DM interaction is changed 
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where J' and D' are the renormalized coupling constants. 
To have a self-similar Hamiltonian, we implement a 7r ro- 
tation around x axis on all sites (of — > —erf , a\ — > —erf)- 
We note to interpret our final results in terms of this 
transformation. The renormalized coupling constants are 
functions of the original ones which are given by the fol- 
lowing equations. 
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, q= \fl + 8D 2 



We will implement this approach in the next sections to 
obtain the phase diagram and entanglement properties 
of the model. 



III. PHASE DIAGRAM 

The RG equations show the scaling of J coupling to 
zero which represents the renormalization of energy scale. 
At zero temperature, a phase transition occurs upon vari- 
ation of the parameters in the Hamiltonian. In the ab- 
sence of DM interaction (D = 0) the ground state of 
the Ising model is the Neel ordered state. However, for 
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D ^ the DM interaction makes a tendency for spins to 
be oriented in the XY plane. A nonzero value of D in- 
creases the fluctuations which destroys the AF ordering 
in z direction at some finite value of D = D c . Simulta- 
neously, the chiral order grows up and will saturate as 
D — > co. The quantum phase transition can be inter- 
preted as the antiferromagnet to saturated chiral (SC) 
order transition at D = D c . The RG flow shows that in 
the AF phase (D < D c = 1), the DM coupling (D) goes 
to zero and in the SC phase D goes to infinity (Fig.©). 
We have probed the AF-SC transition by calculating the 
staggered magnetization Sm (See appendix B) in the z- 
direction as an order parameter (Fig. ([3])), 



Sm — — 
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Sm is zero in the SC phase and has a nonzero value in 
the AF phase. Thus, the staggered magnetization is the 
proper order parameter to represent the AF-SC transi- 
tion. We have plotted Sm versus D in Fig.©. It has 
its maximum value at D = and continuously decreases 
with increase of D to zero at D = 1. Moreover, we have 
calculated the chiral order— (Ch) in the z direction (See 
appendix B) which increases with D and saturates for 
D -» co (Fig©, 
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FIG. 2: (color online) Phase diagram of the Ising model with 
DM interaction. Arrows show the running of coupling con- 
stant under RG iteration. 
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FIG. 3: (color online) Chiral order (filled circles) and Stag 
gered Magnetization (filled squares) versus D. 



The chiral order has a nonzero value in both AF and 
SC phases which can not be a proper order parameter to 
distinguish the quantum phase transition. However, it 
shows that the onset of DM interaction sets up the chiral 
order immediately. A classical picture of the chiral order 
in terms of the spin projection on the xy-plane has been 
plotted in Figgl 

We have also calculated the critical exponents at the 
critical point (D = 1). In this respect, we have obtained 
the dynamical exponent, the exponent of order parame- 
ter and the diverging exponent of the correlation length. 
This corresponds to reaching the critical point from the 
AF phase by approaching D — > 1. The dynamical ex- 
ponent is given by z — [ln( J/ J')d=i]/[^("-s)] — 0.73, 
where ns — 3 is the number of sites in each block. 
The staggered magnetization close to the quntum critical 
point goes to zero like Sm ~ \D — where /3 ~ 1.15 and 

is obtained by /? = [ln(S' M / Sm ^gr? ]1.p=i where 
prime denotes the renormalized quantity. The correlation 
length diverges £ ~ \D — 1| _!/ with exponent v ~ 2.15 
which is expressed by v = \tn(n B )]/ln[ dD d ^ ]\ D= i. The 
detail of this calculation is similar to what has been pre- 
sented in Ref.[li]. 



IV. NUMERICAL RESULTS 

We have implemented the exact diagonalization 
method based on Lanczos algorithm to get the ground 
state properties of the Hamiltonian defined in Eq.([T]). 
The Hamiltonian does not commute with S z = 
(l/2)X^ a fi which imposes to consider the full Hilbcrt 
space for computations. We have considered a periodic 
chain of length N = 8, 12, 16, 20 and 24 in our calcula- 
tions. We have first calculated the ground state energy 
for different sizes. We have observed the size dependence 
for ground state energy is weak which make us to ex- 
trapolate our results to get the ground state energy per 
site (Eq/N) as its thermodynamic limit (N — > co). We 




FIG. 4: (color online) A classical picture of spin orientation in 
the xy plain where the angle between neibouring spins depend 
on the D value (see appendix [Cjl . 
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FIG. 5: (color online) The ground state energy versus D. The 
result of QRG (filled circle) is compared with the extrapolated 
values of exact diagonalization (filled squares). 



FIG. 7: (color online) A closer look at the crossing point 
of Fig. (|6} which shows that different plots do not cross each 
other at a single point. The crossing point of two successive 
size occurs at larger value of D upon increasing size. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 

D 

FIG. 6: (color online) The z-component static structure factor 
at antiferromagnetic ordering vector k = n versus D. A more 
clear picture about the crossing point is plotted in Fig. (|7jl . 



have plotted Eq/N versus D in Fig. ([5]). Our results show 
that the QRG result for E is close to the exact diago- 
nalization one which justifies the correct trend versus D 
although the values have around 10 percent error. It is a 
good evidence that the QRG result is reliable at least to 
get the qualitative picture of the model. 

Since the numerics is done on finite size systems the 
symmetry breaking can not occur in our calculation to 
show the nonzero value of the order parameter. Instead, 
the structure factor shows a divergent bahaviour at or- 
dering momentum by increasing the size of system. The 
structure factor at momentum k is defined by the follow- 
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FIG. 8: (color online) The z-component structure factor for 
different D versus 1/N. Plots for D < 0.8 show diverging 
behaviour as TV — > oo while those for D > 0.8 becomes finite 
in the thermodynamic limit. 



ing relation 

N-l 

S<* a (k) = — J2(<r>i+r)e ikr - (3) 

The ^-component of structure factor (S zz (k)) versus k 
has a sharp peak at k — it for D < 0.8 representing 
the antiferromagnetic order. To justify if the peak corre- 
sponds to the true long range order (LRO) or it is just 
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FIG. 9: (color online) The x-component structure factor ver- 
sus momentum (k) for different D. 
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FIG. 10: (color online) The finite size scaling of S xx {k = tt/2). 
All plot show non-divergent behaviour as N — > oo represent- 
ing no long range order but a local order. 



a local order we have computed the structure factor for 
different size of chains (N). We have plotted in Fig.© 
S zz (k = 7r) versus D for different chain lengths (N). 
The peak height increases for D < 0.8 and decreases 
for D > 0.8. To have a clear picture of these data, we 
have plotted in Fig.® the same data for fixed D value 
versus 1/N. We observe a divergent behaviour for the 
z-component structure factor for D < 0.8 and diminish- 
ing for D > 0.8. This justifies antiferromagnetic LRO for 
D < 0.8. Although the plots for different N in Fig.© 
show to cross each other at a single point D* — 0.76 
the fine tunning data close to this point (Fig. ([7])) rep- 
resent different crossings for two successive N values, 
ft is the manifestation of finite size scaling which ex- 
ist in our numerics. Therefore, the true critical point 
(D c ) which should be the case for N — ► oo is greater 
than D* = 0.76. The investigation which shows the re- 
lation between the fsing model and DM interaction with 
the anisotropic Heisenberg model (XXZ) verifies that the 
critical point should be at D c = 1. ft will be discussed 
in SecED 

To get a picture on the type of ordering in the xy plane 
we have plotted the x structure factor of N = 24 for dif- 
ferent D values versus k in Fig.©. Due to symmetry 
we have S xx {k) = S vy (k), thus, we only present data 
for S xx (k). The x-component structure factor show two 
strong peaks at k = ir/2,3ir/2. It is a justification of 
spiral order in the xy plane. However, this is a local or- 
dering and is not a true LRO. We have also plotted in 
Fig. ([10]) the value of S xx {k = tt/2) versus 1/N for dif- 
ferent D values. All data in Fig. ([i"0|) show diminishing 
behaviour as N — > oo. This justifies that the spiral (chi- 
ral) order which exist in the xy plane for D ^ is not a 
true LRO. 



V. ENTANGLEMENT AND ITS SCALING 
PROPERTY 

In this section we calculate the entanglement of the 
model using the idea of renormalization group2£. As we 
have mentioned previously, a finite size block is treated 
exactly to calculate the physical quantities. The coupling 
constants of a finite size blocks are renormalized via the 
QRG prescription to give the large size behavior. Bipar- 
tite entanglement, i.e the entanglement between some de- 
grees of freedom and the rest of system, is quantified by 
von-Neumann entropy of eigenvalues of the reduced den- 
sity matrix. In our case, we first calculate the entropy of 
the middle site and the remaining sites of a single block 
(see Fig.l). The entanglement is easily calculated, since 
the density matrix is defined by 



Q = IV'oXV'ol 



(4) 



where |"0o) has been introduced in Appendix 
(Eq. (|A2[) '). The results will be the same if we consider 
I^q) to construct the density matrix. 

The density matrix defined in Eq.© is traced over 
sites 1 and 3 to get the reduced density matrix for site 2 
(02) which gives 
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FIG. 11: (color online) Representation ol the evolution of 
entanglement entropy in terms of RG iterations. 



FIG. 12: (color online) First derivative of entanglement en- 
tropy and its manifestation towards divergence as the number 
of RG iterations increases (Fig |ll|) . 



In the spirit of RG, the first iteration of RG repre- 
sents a chain of 3 2 sites which is described effectively 
by three effective sites interacting via the renormalized 
coupling constants. Having this in mind, we understand 
that in the first RG iteration the von-Neumann entropy 
with renormalized coupling constant yields the entangle- 
ment between effective degrees of freedom. The varia- 
tion of entanglement (E) versus D is plotted in FigfTTI 
Different plots show the evolution of E under QRG it- 
erations. In other words, the different iteration of QRG 
show how the entanglement evolves as the size of chain 
is increased. Long wavelength behaviors are captured as 
the RG iterations are increased. In Fig [11] we see that in 
the gapped phase, i.e AF, and the long-wavelength limit 
the entanglement is suppressed while in the SC phase the 
entanglement gets maximum value due to the DM inter- 
action in the XY plane that induces a state with strong 
quantum correlation. Such a behavior has also be seen 
in the XX Z model^l. 

A common feature of the second order phase transi- 
tions is the appearance of nonanalytic behavior in some 
physical quantities or their derivatives as the critical 
point is crossed^. It is also accompanied by a scal- 
ing behavior since the correlation length diverges and 
there is no characteristic length scale in the system at 
the critical point. Entanglement as a direct measure of 
quantum correlations indicates the critical behavior such 
as diverging of its derivative as the phase transition is 
crossed^. It has been verified that the entanglement in 
the vicinity of critical point of Ising model in transverse 
field (ITF) and XX model in transverse field shows a 
scaling behavior. Investigating the nonanaliticity, e.g a 
divergence, and finite size scaling provides excellent esti- 
mates for the quantum critical point. A precise connec- 
tion between the entanglement in quantum information 
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FIG. 13: (color online) The scaling behavior of D m in terms 
of system size (N) where D m is the position of minimum in 
Figrj2l 



theory and the critical phenomena in condensed matter 
physics has been established^, where the scaling proper- 
ties of the entanglement in spin chain systems, both near 
and at a quantum critical point have been investigated. 
The first derivative of entanglement let us to get more in- 
sight on the qualitative variation of the ground state as 
the critical point is touched. To this end we have calcu- 
lated the first derivative of entanglement which has been 
depicted in Fig |121 Such a computation determines the 
scaling law of entanglement in one-dimensional spin sys- 
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FIG. 14: (color online) The logarithm of the absolute value of 
minimum, \n(dE/dD |o m ), versus the logarithm of chain size, 
ln(iV), which is linear and shows a scaling behavior. Each 
point corresponds to the minimum value of a single plot of 
FigUll 



FIG. 15: (color online) The finite size scaling is performed 
via the RG treatment for the power-law scaling. Each curve 
corresponds to a definite size of the system, i.e N = 3™ . 
The exponent 9 is ascribed to the correlation length critical 
exponent v via 9 = 1/v 



terns, while explicitly uncovering an accurate correspon- 
dence with the critical properties of the model. As the 
size of system becomes large through RG iterations, the 
derivative of entanglement tends to diverge close to the 
critical point. All plots in FigQ2]with respect to the crit- 
ical point have an asymmetrical shape. Each plot reveals 
a minimum in the gapped phase, i.e AF for < D < 1, 
the minimum becomes more pronounced close to the crit- 
ical point, D = 1. It manifests that the ground state of 
the gapped phase of the model undergoes a strong qual- 
itative change when approaching the quantum critical 
point while the corresponding change in the SC phase is 
rather small. A similar situation has also been observed 
in the XX Z model^i. This behavior is comparable with 
results on ITF model where the system in both sides of 
the critical point is gapfull, so the derivative of the en- 
tanglement tends to diverge symmetrically^. 

More information can be obtained when the minimum 
values of each plot and their positions are analyzed. The 
position of the minimum (D m ) of ^ tends towards the 
critical point like D m = D c — N~- aA6 which has been 
plotted in Fig[l3l Moreover, we have derived the scaling 
behavior of y = |^|c m versus N. This has been plotted 
in Fig ll4l which shows a linear behavior of ln(y) versus 
ln(N). The exponent for this behavior is | 4§|.d„, ~ 
A 46 . This results justify that the RG implementation 
of entanglement truly capture the critical behavior of the 
model at D = 1. It should be emphasized this exponent 
is directly related to the correlation length exponent, v, 
close to the critical point. It has been shown in Ref.|27| 
that | % \ Dc ~ N xjv and D m = D C + N~ x ' v . 

To study the scaling behavior of the entanglement en- 



tropy around the critical point, we perform finite scaling 
analysis. Since the minimum value of derivative of en- 
tanglement entropy scales power-law. According to the 
scaling ansatz, the rescaled derivative of entanglement 
entropy around its minimum value, D m , is just a func- 
tion of rescaled driving parameter like: 

dE _ dE I 

dD N f lDm — F[N e (D — D m )] 

where, F(x) is a universal function that does not deponed 
on the system size, and the exponent 9 is just the inverse 
of the critical exponent v, i.e 6 = The manifestation 
of the finite size scaling is shown in FigJT5] It is clear 
that the different curves which are resemblance of various 
system sizes collapse to a single universal curve. It must 
be noticed that the n—th RG iteration describes a system 
with 3 ,l+1 sites which is effectively represented by a three 
site model through the RG treatment. 



VI. SUMMARY AND CONCLUSIONS 

We have applied the quantum RG approximation to 
obtain the phase diagram, staggered magnetization, chi- 
ral order and the entanglement properties of Ising model 
with DM interaction. Tunning the DM interaction dic- 
tated the system to fall into different phases, i.e antifcr- 
romagnetic phase with nonzero staggered magnetization 
(as order parameter) and chiral one with vanishing order 
parameter. The critical point of the phase transition is at 
D c = 1 where the quantum fluctuations have dominant 
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effect which arises from the DM interaction and eventu- 
ally destroy the order in the AF phase. Although the 
DM interaction drives the spins to leave their ordering 
in the z-direction, i.e staggered magnetization, a chiral 
saturated phase has been arisen, Figf2j 

The numerical exact diagonalization which has been 
done on finite sizes (N — 12, 16, 18, 20) justifies the re- 
sults of QRG. We have obtained a fairly well agreement 
in the ground state energy of QRG approach with exact 
diagonalization method. The divergence of ^-component 
structure factor at momentum k = n when N — > oo 
is the signature of antiferromagnetic long range order 
for D < D c . The value of critical point which is read 
from exact diagonalization is £>^ umerlc ~= 0.8 which is 
20 percent different from what QRG gives, D^ RG = 1. 
However, we claim that the QRG result for D c should 
be more reliable as will be discussed below in connection 
with anisotropic Hciscnbcrg (XXZ) model. Moreover, the 
size dependence of the critical point (D c ) is strong which 
demands the numerics on larger sizes and also a scaling 
analysis. 

Besides, the entanglement entropy of the model at dif- 
ferent RG iterations was analyzed. As the long wave- 
length behavior of the model is reached via the increasing 
of the RG iterations, the entanglement entropy develops 
two distinct behavior proportional to two different exist- 
ing phases of the model. However nonanalytic behavior 
close to the critical point of the model manifests itself 
via the analysis of the first derivative of the entangle- 
ment entropy. The divergence of the first derivative of 
the entanglement entropy becomes more pronounced as 
long as the size of the system becomes large in RG treat- 
ment. Critical point is touched by an exponent which ap- 
pears as the inverse of the critical exponent which shows 
the divergence of the correlation length. Moreover, it 
is found that as the critical point is touched from the 
gapped phase a drastic change in the ground state oc- 
curs which manifests itself in the evolution of the deriva- 
tive of entanglement (see Fig fT2")l . Such variation in the 
ground state structure also appears in the XXZ model. 
Finally finite size scaling reveals the critical properties 
of the model is mirrored via the nonanalytic behavior of 
the entanglement. 

The one dimensional Ising model with DM interaction 
(Eq.([T])) is mapped to the XXZ chain via a nonlocal 
canonical transformatio n 36 ! 37 , 
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j'-i 



tj = to tan 1 (D), 



of = e- iU afe iU , a* = a z , 
H = e- lU He tU , 



(7) 



which gives 

H ~ £ (dfSf ¥1 + a\a\ +1 + ih^t+x)- (8) 



This transformation tells that the system is in the 



spin-fluid phase for D > 1 in terms of transformed 
spins. While, the model represents the antiferromagnetic 
(Neel) phase for D < 1. The spin- fluid phase (D > 1) 
of the transformed Hamiltonian (H) can be character- 
ized by a string order parameter—. At the same time the 
saturated chiral phase (D > 1) of the original Hamilto- 
nian (Eq.(fT])) is represented by the chiral order. It can 
be concluded that the chiral order with nonlocal spins 
is similar to the string order parameter. Therefore, it is 
suggested that the XXZ model has chiral order which is 
constructed with nonlocal spins. It will be instructive to 
calculate the chiral order with nonlocal spins for XXZ 
model as a hidden order of the spin-fluid phase. 
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APPENDIX A: THE BLOCK HAMILTONIAN OF 
THREE SITES, ITS EIGENVECTORS AND 
EIGENVALUES 

We have considered the three-site block (Fig.(fT])) with 
the following Hamiltonian 



hf 



J r 



<A + °h a li - ■ ( A1 ) 



The inter-block (H BB ) and intra-block {H B ) Hamil- 
tonian for the three sites decomposition are 



H 1 



N/3 

E 

i=l 



~ z ~ z i ~ z ~ z 
Y CT 1,I (J 2J + (J 2.I (J 3J 



J 

4 ^ 

=1 



H 



BB 



N/3 

-Y 

i=i 



where a"j refers to the a-component of the Pauli matrix 
at site j of the block labeled by /. The exact treatment of 
this Hamiltonian leads to four distinct eigenvalues which 
are doubly degenerate. The ground, first, second and 
third excited state energies have the following expressions 
in terms of the coupling constants. 



9 



\i'2 



('3 



m = 


1 




V * v 


+ 9) 


Wo) — 


1 




V 2 9(l 

V * \ 


+ ?) 


eo = 






l^i) = 


1 






-1) 


Wi) = 


i 




^29(9 


-1) 


ei = 


-£(i- 


9), 



[2£>| 1TT)+*(1 + ?)I TIT) 
[2£>| UT>+*(i + ?)l TTT> 

[2D\ lTT)-*(«-i)l TTT) 
[2D\ Ut) -*(«-i)| TTT) 



2£>| TTT)] 
2£>| TTT)] 



function at each iteration of RG is connected to its value 
after an RG iteration. This will be continued to reach a 
controllable fixed point where we can obtain the value of 
the correlation function. The staggered magnetization in 
a direction can be written 



S 



1 N 



(-i)* 



(Bl) 



2D\ TT1)]) wnere a f is the Pauli matrix in the ith site and \0) is the 
ground state of chain. The ground state of the renormal- 
ized chain is related to the ground state of the original 

2-D| TU)J, one by the transformation, Pq\0') = \0). 



N 



^=(1 TTT) + I TTT)) 



W) = ^(l TTT) + I TTT)), 



N 



= 0, 



('2 

m = I TTT) , |V>3> = I TTT), 



J 
2' 



This leads to the staggered configuration in the renormal- 
ized chain. The staggered magnetization in z direction is 
(A2jbbtained 



where q is q = vT + 8D 2 . 

I j) and I |) are the eigenstates of o z . The projection 
operator is 

Po = | Wo| + I Wol- 

The Pauli matrices in the effective Hilbert space have the 
following transformations 



P 1 rr x P 1 
pi x pi _ 

pI n V pi 
M) a 2,l — 



P 1 rr z P 1 



P 1 rr z P 1 



2D 



IV pi x pi 



4L> 2 



2D 



/.v 



cr 

4I? 2 

9(9 + 1)" 
i + g ,z 

~2q~° 1 



P 1 ^ P 1 



9(9 + 1) 

2D 



1 , 



M) cr 3,7- r 



7 

2D 



1 i 



p'n- 2 P j - --rr' 
M) a 2,I r Q — „" I > 
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APPENDIX B: ORDER PARAMETER AND 
CHIRAL ORDER 



N 



(-1)* 



^10) 



1 1 

6 



W/3 



3 7=1 

/+!/ 



(0'|P J+1 (-<x + i + ^ + i 

JV/3 



)^E<° 



-<7)^0|0') 

^,/ + i)^o +1 |0') 

7° 



3 ? 



where is the staggered magnetization at the nth step 



of QRG and 7W is defined by 7 



9)/9- 



This process will be iterated many times by replac- 
ing 7^°) with 7W. The expression for is similar to 
7C ) where the coupling constants should be replaced by 
the renormalized ones at the corresponding RG iteration 
(n). The result of this calculation has been presented in 
Fig.©. 



1. Staggered magnetization 



Generally, any correlation function can be calculated 
in the QRG scheme. In this approach, the correlation 



2. Chiral Order 



10 



c h = 



^Ei<oi(^r+i-^i)io) 

i=l 
1 1 N/3 

Y^Z £ [<°'I P °KX /+1 ~ <z</+i)^o|0') + (0'|itf ((of,,**/ - o^of,,) + (of,,**, - o^at tI ))Pi |0') 



1 1 

Tv 

3 /=1 

1 32L> 3 
12? 2 (l + g) 



3 J — l 



1 32D 3 

? 2 (i + g) 



T° = (^) 2 . 



At the last step we use the following transformation, 
u « ' u i u i ■ 



projection on the z-axis and have the helical structure 
(see Fig|4]) in the xy plain. In this region the minimum 
classical energy is 



APPENDIX C: CLASSICAL APPROXIMATION 

In the classical approximation the spins are considered 
as classical vectors 

erf = cos{ip) sva.9 , of = sm(iip) sin (9 , erf = cos#, 

where ip is the azimuthal angle measured from the x- 
axis in the xy-plane and 9 is the polar angle measured 
down from the z-axis. The classical energy per site for 
IDM Hamiltonian (Eq.©) is: 



E c i 
N 



J 



(cos 2 9 + D sin 2 9 sin ip). 



The minimization of classical energy with respect to the 
angles ip and 9 shows that there are two different regions. 
(I) D > 1, the minimum of energy is obtained by 9 = 5 
and <p = arcsin(-^j) which show the spins do not have a 



N 



J 

T 



(Cl) 



(II) D < 1, the energy is minimized by (p = § an< A 
arbitrary 9 which corresponds to the configuration with 
nonzero value of spins projection on z-axis and helical 
structure of spins projection in the xy-plain where the 
angle between spins are ^. In this region the minimum 
classical energy is 



El = i(cos 2 9 + Dsin 2 9). 
N 4 v ' 



(C2) 



One can see from Eq. (|Cljl and Eq. (|C2j) that the tran- 
sition between phase (I) and (II) takes place at D = 1. 
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